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We show that the propagation of spin polarized carriers may be dramatically affected by inhomogeneous 
electric fields. Surprisingly, the spin diffusion length is found to strongly depend on the sign and magnitude 
of electric field gradients, in addition to the previously reported dependence on the sign and magnitude of 
the electric field [Yu & Flatte, Phys. Rev. B 66, 201202(R) 2002; ibid. 66, 235302 (2002)]. This indicates 
that purely electrical effects may play a crucial role in spin polarized injection, transport and detection in 
semiconductor spintronics. A generalized drift- diffusion equation that describes our findings is derived and 
verified by numerical calculations using the Boltzmann transport equation. 
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I. INTRODUCTION 

Semiconductor spintronics has recently received a lot 
of attention following the successful demonstration of 
spin injection from a diluted magnetic semiconductor 
(DMS)^ to a nonmagnetic semiconductor (NMS).^ This 
enables the realization of all-semiconductor spintronics, 
which potentially can lead to the development of multi- 
functional devices based on optical and transport pro- 
cesses that make use of both the charge and spin de- 
grees of freedom. However, in order to realize spintronic 
devices, several prerequisites need to be generally satis- 
fied and achieved: i) efficient creation and injection of 
spin polarized carriers, ii) transport of the spin-polarized 
electrons from one location in the device to another, iii) 
successful manipulation of the spin current signal, iv) ef- 
ficient detection of the spin-polarized electrons. 

While the electron spin is of central importance in the 
understanding and design of semiconductor spintronics, 
the charge of the electron and its influence on spin trans- 
port cannot be neglected. Charge carrier interactions, 
drift and diffusion are expected in any realistic system 
due to externally applied, or intrinsic, built-in electric 
fields, doping concentration variations and material de- 
sign. It has been previously shown by Yu and Flatte (YF) 
that the spin diffusion length depends on the magnitude 
and sign of the electric field J- Such a dependence was re- 
cently experimentally observed^ It has also been shown 
experimentally that spin-polarized transport and injec- 
tion is very sensitive to an applied bias volt age f^i^i^i^ 
as well as different structural parameters such as doping 
concentrations^ and layer thicknessesj^ These considera- 
tions pose the fundamental question: What is the influ- 
ence of purely electrical effects on spin injection, trans- 
port and detection in semiconductor spintronics? 

Here we show that the the spin diffusion length may 
have a dramatic dependence on electric field gradients, an 
effect that can be significantly stronger than the electric 
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field effects predicted by YF.~ In essence, our theoreti- 
cal and numerical study predicts that even small devia- 
tions from charge neutrality can give rise to significant 
enhancement or suppression of spin polarization. 

Previously, it has been assumed that the effects of in- 
homogeneous fields, such as, e.g., occurring at material 
or doping concentration interfaces, should be relatively 
small, since the charge screening length typically is much 
shorter than the spin diffusion length. We show that the 
effective spin diffusion length in fact can be comparable 
to the screening length in the presence of inhomogeneous 
electric flelds. In particular, we study the propagation of 
spin polarized electrons across interfaces between regions 
with different doping concentrations. The spin polariza- 
tion of electrons moving from a region with lower to a 
region with higher doping concentration is found to be 
strongly suppressed due to the electric field gradient ef- 
fect. Conversely, electric field gradients may enhance the 
spin diffusion length, depending on their sign and mag- 
nitude. 

Our results show that purely electrical effects may have 
a profound impact on spin injection, spin transport, and 
spin detection, three key ingredients for the realization of 
semiconductor spintronics. A generalized drift-diffusion 
equation which is able to describe nonequilibrium spin 
density propagation in the presence of inhomogeneous 
electric fields is derived and verified with numerical calcu- 
lations using the Boltzmann transport equation (BTE). 
We note that recent theoretical studies by, e.g., Wang 
and Wu^ and Sherman^ have also predicted very inter- 
esting effects related to electric-field induced spin relax- 
ation mechanisms due to the Rashba effect. 

In the following, we will first describe our self- 
consistent numerical approach (Section |TT| , followed by 
numerical results from which the strong influence of elec- 
tric fields on spin transport will become evident (Sec- 
tion IIII[) . Subsequently, in Section IIVI we derive spin 
drift-diffusion equations which are then used to discuss 
and understand the numerically observed characteristics 
(Section |V|. Finally, our conclusions are presented in 
Section IVTl 
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II. SELF-CONSISTENT BOLTZMANN-POISSON 
MODEL 

In order to fully understand spin transport properties 
in semiconductor structures, and in particular the influ- 
ence of interfaces and inhomogeneous electric fields, we 
take into account nonequilibrium transport processes for 
both the charge and spin degrees of freedom in a self- 
consistent way. In our approach, the transport of spin- 
polarized electrons is described by two BTE equations 
according to 
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where E is an inhomogeneous electric field, fm) is the 
electron distribution for the spin- up (down) electrons. 
Tin is the momentum relaxation time, and (l/''"iT) 
is the rate at which spin-up (spin-down) electrons scat- 
ter to spin-down (spin-up) electrons. The first term on 
the right-hand side of eq. ^ describes the relaxation 
of each nonequilibrium spin distribution to a local equi- 
librium (spin-dependent) electron distribution function, 
/|^-|, which we choose as non-degenerate according to 
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where T is the lattice temperature. The last term in 
eq. ([T]) describes the relaxation of the spin polarization. 
From the distribution function we calculate the local 
spin density according to 



(3) 



T(i)W = / /T(i)(i"'V)dv 



and define two quantities, the spin density imbalance 

= "T - "i ' (4) 

where and are the densities of spin-up and spin- 
down polarized electrons, and the spin density polariza- 
tion 



Pn = 



n 



(5) 



where n = + is the total charge density. 

Inhomogeneous charge distributions and electric fields 
are taken into account by solving the Poisson equation 



-V • E 



Nd 



eso 



(6) 



where e is the dielectric constant, Nd is the donor con- 
centration profile, and (f> and E are the spatially depen- 
dent potential and electric field profiles. 

Our theoretical model thus goes beyond drift-diffusion, 
and is capable of describing charge and spin transport 
through strongly inhomogeneous (with respect to spin 
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FIG. 1: (color online) Schematics of the model structure we 
consider. Full spin polarization, Pri = 1 is assumed for x < Xp. 



and charge densities) semiconductor systems, as well as 
nonequilibrium effects. We assume in the following that 
an electric field is applied along the x direction and con- 
sider the corresponding one-dimensional transport prob- 
lem. 



The nonlinear, coupled equations (Il|2|3l6p need to be 
solved self-consistently for a given applied bias voltage, 
doping concentration profile, nonequilibrium spin distri- 
butions, and relaxation times. We use a numerical ap- 
proach based on finite difference and relaxation methods, 
that wc originally developed for the study of nonequi- 
librium effects in charge transport through ultrasmall, 
inhomogeneous semiconductor channels .-iiiJ^ As bound- 
ary conditions, we adopt the following scheme: For the 
potential, the values at the system boundaries are fixed 
to (j){xi) = Vb and (j){xr) = (^,7- denote the left and 
right boundary of the sample, respectively), correspond- 
ing to an externally applied voltage Vt- The electron 
charge density is allowed to fluctuate in the system sub- 
ject to the condition of global charge neutrality, which 
is enforced between each successive iteration in the self- 
consistent Poisson-Boltzmann loop. The spin density at 
the DMS boundary and in the DMS is determined by the 
degree of polarization P„ = (n-j —n^)/ (n^ +ni), for which 
the density is defined according to »t-|(x) = n/2(l ± P). 
In the calculations, the size of the contacts has to be 
large enough, such that the electric field deep inside the 
contacts is constant and low. This allows us to use the 
analytical, linear response solution to the BTEs ^ 

fKl){^Lr,v) = f'^i^^^{xi^r,v)[l - vE{xi^r)Tm/kBT] , (7) 

as boundary conditions at xi^r, where we use the local 
equilibrium distribution f'^^^^(xi^r,v) and local electric 
field, E{xi^r), obtained from the previous numerical so- 
lution to the Poisson-Boltzmann iterative loop. More 
details of the numerical procedure for a similar problem 
are given in Ref. ,13) . 

The model structure that we study is described in Fig. 
[H As seen in the figure, the system is inhomogeneous 
with respect to both charge and spin degrees of freedom. 
Spin-polarized electrons, with P„ = 1, are injected from 
a diluted magnetic semiconductor(DMS)-like portion,—— 
defined for x < Xp, into a nonmagnetic semiconductor 
(NMS) part. The spin-polarized electrons subsequently 
relax due to spin-flip scattering for x > Xp, at the rate 
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I/ts/ (where t^i = r^i = 2ts/), to the asymptotic unpo- 
larized value P„ = for x ^ Xp. The doping concentra- 
tion is defined as Nd = Nl and Nd — Nu for x < Xd and 
2; > Xd, respectively. The following material and system 
parameters have been used in our study: T = 300 K, 
m* = 0.067mo, r,n — O.f ps (typical for GaAs at room 
temperature), Tg/ = 1 ns, Nl = 10^^ m^'^, Vb = —0.3 V. 
The total system size is 5 /xm, i.e., xi^r = t2.5 /im. 



III. NUMERICAL RESULTS: ENHANCEMENT 
AND SUPPRESSION OF SPIN POLARIZATION 
AT DOPING INTERFACES 

In the following, we present numerical results for a 
spin and charge inhomogeneous semiconductor structure 
and highlight the characteristics of electric field effects 
on the spin polarization. In Fig. [2] we study the prop- 
agation of spin-polarized electrons across two spin- and 
charge-inhomogeneous interfaces, separated by L — 0.2 
/im (grey mid-region), as a function of different doping 
concentrations Nu. In Fig. [2ja) we show the spin den- 
sity imbalance, S-^i = — (solid lines), and the total 
charge density, n (dashed lines), for Nji = rN^, where 
r = 1, 2, 5, 10, 20 {Nl = lO^i m^^. The lowest solid and 
dashed curves correspond to d-\i and n of a charge ho- 
mogeneous sample, Nr = Nl, in which the electric field 
is accordingly constant, as shown by the thick, solid line 
in Fig. [2jd). In contrast to the constant total charge 
density (bottom dashed line), for the charge homoge- 
neous case, the spin density imbalance, S-^i (bottom solid 
line), decreases monotonically for x > Xp, which coin- 
cides with the position of the DMS/NMS interface and 
the onset of the spin flip scattering rate 1/ts/. The de- 
cay is exponential and can be understood in terms of the 
drift-diffusion description in the YF modeli^ Within this 
model, the spin density imbalance decay can be described 
by S^i exp[-a;/L£,([j)] where 
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and where, ^^(m) are electric-field dependent spin dif- 
fusion lengths, describing spin propagation antiparallcl 
(parallel) to the electric field. In eq. ^ Lg = ^Dtsj 
is an intrinsic spin-diffusion length in the absence of an 
electric field and D = ksTTm/m* , as obtained from the 
Einstein relation. From eq. ^ it follows that the spin- 
diffusion length is enhanced in the direction anti-parallel 
to an applied electric field and suppressed in the direc- 
tion parallel to the field.'' Our numerical results agree 
with these considerations for the homogeneous system. 



For increasing Njj, the total charge density, n, nat- 
urally increases monotonically for x > x^, reaching 
n « Nn on the order of a screening length to the right 
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FIG. 2: (color online) Spin density imbalance Sn (solid lines) 
and total charge density n (dashed lines) for a spin and charge 
inhomogeneous structure as depicted in Fig. [1] for different 
doping concentrations (a) Nr = 1,2,5,10,20 • Nl, and (b) 
Nr — 0.05, 0.2, 0.5, 0.8, l-A'i, (increasing from bottom to top), 
(c) Spin density polarization P„, and (d) electric field distri- 
bution for the same structure. The structure parameters are 
Nl = 10^^ m~^ L = 0.2 /im, r^/ = 0.5 ns, = 0.1 ps, 
Vb = -0.3 V. 
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of the doping interface. The spin density imbalance Sn, 
however, has a more comphcated behavior. At the in- 
terface between the two doping regions, defined at Xd, 
a sharp increase of is observed, which subsequently 
peaks below the maximum value of n, followed by a 
monotonic decay. Similar spin accumulation was recently 
reported by Pershin and Privman.— We note and empha- 
size, however, that the spin accumulation peak is signif- 
icantly lower than the maximum value of n around the 
doping interface, e.g., for the Nji = 20-10^^ m~^ sample, 
Sii reaches only « 50% of the corresponding value in n. 

For Nji < Nl, when electrons are injected from a high 
to a low doping region, the situation is reversed. In Fig. 
[2ljb), the dashed and solid lines correspond to n and 
Sii, calculated for different values of Nn = rN^, where 
1 > r > 0.05. For decreasing Nn, S-^ (and naturally n) 
decreases as well. 

The observed increase (decrease) of S^i with increas- 
ing (decreasing) Nn can at first order be understood in 
terms of the large difference between the intrinsic spin 
relaxation length, L^, and the Debye screening length, 
L]jB- The screening length is generally shorter than the 
spin relaxation length. For example, for the range of 
doping concentrations considered here, we find 0.03 < 
Ldb < 0.6 /im, whereas Ls ~ 1.8 /im. The spin density 
imbalance can be written in terms ofSn = n— 2ri|. Since 
Ls > Ldb, the total charge density, n, will increase (de- 
crease) faster than n<^ with increasing (decreasing) Nr. 
Hence, one would expect the spatial dependence of n to 
dominate (5f j within a screening length of the doping in- 
terface. This is clearly seen by comparison between S^i 
(solid lines) and n (dashed lines) in Figs. [2l^a,b). 

For Xd, ^ix is s^^^ to decrease exponentially. Fur- 
thermore, the decay rate is larger for increasing doping 
concentrations Nji. One can understand this spatial de- 
pendence from previous arguments and the model of YF. 
For X ':S> Xd the total charge density n and the electric 
field is constant (and thus we can apply the YF model) 
and the spin-flip scattering drives the noncquilibrium 
spin density distribution exponentially toward equilib- 
rium, 5^1 = on a length scale given by eq. ([8]). Since 
an increase in Nn gives rise to a decrease of the electric 
field in that region [see Fig. Hfd) for the electric profiles] 
the spin diffusion length correspondingly decreases con- 
sistent with the predictions of eq. giving rise to a 
faster decay of the spin density imbalance. 

The description provided by eq. ([8]), however, fails to 
describe the magnitude, as well as the spatial dependence 
of around the doping interface. For example, consid- 
ering first the curves with Nn > Nl, for Nn = 2Nl, 
5-\i follows the spatial profile of n rather closely around 
the interface region. Surprisingly, with increasing Njf , 
the maximum value reached by (5f | becomes increasingly 
smaller than N^, and the spin density imbalance thus 
gets strongly suppressed. This is an unexpected find- 
ing, and in contradiction the YF model and with eq. ([8]), 
since, according to this equation, a strong built-in nega- 
tiv field, such as shown in Fig. lUd), around the doping 



interface should enhance the spin diffusion length signif- 
icantly and hence, should yield (5|| w n until the electric 
field drops again for x >> Xd- 

This discrepancy is also well illustrated in the spin den- 
sity polarization, P„, shown in Fig. [DJc). Although, as 
naively expected, P„ = Sn/n drops sharply around the 
interface due to the sharp rise in n, the dependence on 
Nj^ is in contradiction with eq. ^ due to the following: 
An increase in iV^, yields a sharp increase of around 
the interface [see Fig. [2Kd)], and should therefore lead to 
a significant increase of the " downstream" spin diffusion 
length, Ljj, according to eq. (H]). Thus, an increase in 
Nn should be beneficial for S-^i and P„, in contradiction 
to our findings. For < Nl, a similar unexpected 
dependence on Nj^ is observed. 

Clearly, one needs to go beyond the YF model (which 
is valid for a system with local charge neutrality) in order 
to understand spin polarized transport in the presence of 
inhomogeneous electric fields . As we will discuss below, 
the sign and magnitude of the electric field and electric 
field gradients around the interface are crucial for deter- 
mining the magnitude and spatial dependence of the spin 
density imbalance and polarization around the doping in- 
terface. 



IV. SPIN DRIFT-DIFFUSION EQUATIONS 

A qualitative understanding of the unexpected spa- 
tial dependence of P„ and Sn, can be obtained by re- 
examining the spin drift-diffusion equations. Previous 
works have mainly focused on the charge homogeneous 
casCi ^'^^i^^'^^ This is a natural assumption since the 
screening length is expected to be much shorter than the 
intrinsic spin relaxation length, leading one to assume 
quasi-local charge neutrality. However, as we will show 
below, corrections to the spin drift-diffusion equations of 
YF- yield that even quasi-local deviations from charge 
neutrality can have a very strong impact on the spin po- 
larization. 

In the following, we perform the same steps as YF— 
for the derivation of the spin drift-diffusion equations. 
The current for spin-up and spin-down electrons can be 
written as 

hil) = f^T(i)E + eD|(x)Vn|(|) , (9) 

where ^|(x) and cr|(x) are the diffusion constants and 
conductivities for the spin-up(down) species. Here, the 
conductivities refer to crm) = enxQ-j/xm), in terms of 
the mobilities Hid) ^^'-^ densities n^d) of spin-up(down) 
electrons. The continuity equations for the two species 
require that 

r.. = \ ^-V-JX(X), (10) 

Ot ni(xx) Txx(Ti) e 

where T||(t||) are the spin-flip scattering rates. In 
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FIG. 3: (color online) Total charge density (dashed line) and 
spin density imbalance (solid line) profiles calculated for Nr — 
20Nl- Additional curves for —0.1 < a; < 0.1 /im correspond 
to 5^1 = n{x) exp[— (a; + 0.1)/LD] for Ld — 0.16, 0.3, 1, 36 /im, 
as described in the text. Inset shows n and (5|| for a larger 
region. 



steady-state eqs. (|9ll0p result in 

VCT|(x)-E+(7|(|)V-E+ei:i|(x)V2n|(|) 

^Ti(iT) nT(Ti) 
(11) 

For a NMS = ^4, and ^ D^. Multiplying eq. 
(fTTj) with a I (cTj^ ) and substracting them from each other 
we arrive to the following expression for the spin density 
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where we have used the Einstein relation fi — eD /ksT^ Tg 



is the spin relaxation time r^, ^ = 
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and Lg = yDrl is the intrinsic spin diffusion length. 



Equation ([T2|) is a drift-diffusion equation for the spin 
density imbalance and contains terms depending on 
both the electric field, and the electric field gradient. We 
note that in fact eq. p2|) resembles the corresponding 
equations in Refs. 3, but with an additional term pro- 
portional to V • E. Similarly to Refs. [Sj we can use the 
roots of the characteristic equation for eq. ([T^ to define 
up- and down-stream spin diffusion lengths from eq. (jl2p 
leading up to the following equation 



L' 



D{U) 



-( + ) 



\eE\ 
2kBT 



(2fcsr) 



1 



eVE 



(13) 

Notice that eq. (fT51) is only defined "locally" , over a re- 
gion where V • E can be considered constant, and where 



an average value of the electric field is used. Note also 
that the above equations are also valid for ferromagnetic 
semiconductors provided the mobility and diffusion con- 
stant are replaced with the ones corresponding to the 
minority-spin speciesi^ii^ 

A comparison between eqs. and ([5]) shows that 
eq. (|13p contains an additional term —eVE/ksT in the 
square root expression. This term is a correction to the 
YF model which describes the important role that in- 
homogeneous electric fields can play on spin polarized 
transport. 



V. DISCUSSION 

From eq. it is thus seen that the drift-diffusion 
length depends strongly on the gradient of the elec- 
tric field; Consequently, the spin density imbalance and 
polarization decay lengths are enhanced or suppressed 
based not only on the sign of the electric field alone 
On the contrary, depending on the profile of the electric 
field, the gradient term in eq. (jl3p can either enhance or 
suppress the spin-polarization propagation significantly, 
as seen explicitly in our calculated results in Fig. [2I^a-c). 

To illustrate the dramatic impact of the electric field 
gradient term on the spin polarization diffusion, we cal- 
culate i^((7) for the central region of the structure, 
Xp < X < Xd- Using a linear approximation of the 
electric field within this region, for the structure with 
TVfi = 20 • 10^1 m-3 [bottom, sohd fine in Fig. ^c)], the 
intrinsic spin-diffusion length Lg sa 1.8 /Ltm, and an elec- 
tric field value taken at mid-channel, E = —2.9 kV/cm, 
eq. yields L'^ « 0.16 ^m and L'jj « 0.06 /xm. In com- 
parison, an evaluation of eq. ([5]) using the same value of 
the electric field, yields Ld k, 36 /im, Lu k, 0.09 /im. 
Hence, the electric field gradient in this case decreases 
the "local" (around the interface) spin diffusion length, 
along the direction of (charge) transport, by two orders 
of magnitude! The injected spin density polarization is 
thus destroyed by the inhomogeneous electric field at the 
doping interface. 

To verify the validity of eq. (|13p . we compare exp[— (x— 
using L'j-, — 0.16 /im as before, with the results 
(for P„) from our Boltzmann-Poisson numerical calcula- 
tion, within Xp < X < Xd- The results, plotted in Fig. 

(thick dashed line for the highest Nj^, value), are 
found to agree very well with the numerically calculated 
values. 

We further do a similar analysis for the spin den- 
sity imbalance, <5||. Here, we compare the numerically 
calculated values for 5^i, and n for the same doping 
concentration profile with = 20Nl, with S'-^^ = 
n{x) eyiY>[—{x — xq) / L'jj], using the numerically calculated 
values of n{x), and the estimates of Ljj and L'jj above. 
The results are shown in Fig. [3]for — 0.16, 0.3, 1 and 
36 /zm, the first and latter values corresponding to the 
previously estimated values using eqs. [T3] and [5] The 
numerically calculated curve of Sn is best fitted to an 
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effective spin diffusion length of L'jj = 0.3 /im, which 
is in good agreement with the value obtained from eq. 
(fT3| . any discrepancies being due to the approximations 
in the choice of VE and E in eq. p^ . Nevertheless, it is 
clear that the spin diffusion length around the interface 
is much smaller than the intrinsic one, Ls, and the cor- 
responding L/j, obtained within the YF model, which is 
only valid in the (locally) charge neutral case. 

Reversing the sign of the electric field gradient, the 
spin density polarization P„ is in fact enhanced in com- 
parison to the homogeneous case. This is illustrated by 
the dashed-dotted curves in Fig. [5]^c), which correspond 
to iVfl < Nl. Since the electric field gradients are not 
as large as for the previously studied case [see Fig.jSlJd)], 
the corresponding increase in P„ is more modest. 

We further clarify the effects of electric field gradients 
on the relaxation of nonequilibrium spin distributions. 
Figure m shows the results for (a), P„ (b), and E{x) 
(c) calculated for Nr = IQNl = 10^^ m'^ and differ- 
ent channel lengths, < L < 3 /xm. We have offset 
all curves in space such that Xp, i.e., the position of the 
DMS/NMS interface, coincides. For i = 0, 5-^i shows 
a sharp increase around the doping interface, followed 
by a monotonic, exponential decrease, whereas the spin 
density polarization, P„, shows an exponential decrease 
only, with decay length similar to (5||, but with values be- 
low the corresponding values for the homogeneous system 
(thick solid line) . This is in agreement with our previous 
discussion and agrees with the fact that the electric field 
[see Fig.|4l[c)] in the right-hand side region is significantly 
lower than for the homogeneous case [thick solid line in 
Fig. He)]. 



For increasing channel lengths the spin imbalance de- 
creases significantly. In particular, 6'\i for the inhomo- 
geneous structure is smaller than the corresponding ho- 
mogeneous one (thick solid line) for x < Xd and a dip- 
like feature is formed for large L prior to the interface 
between the two doping regions, where is seen to in- 
crease and reach a peak. Interestingly, this also occurs 
for X < Xd — LuB for the structures with L 3> L^b, 
i.e., in a region with very small electric field gradients 
[see Fig. IH^c)]. This illustrates the fact that, although 
the electric field gradients in this region are very small 
[virtually unobservable on the scale shown in Fig. HKc)], 
they are large enough to yield significant influence on the 
spin density imbalance. In the spin density polarization, 
Pn, the inhomogeneous electric fields manifest as nonex- 
ponential spatial dependence, which further accentuates 
the importance of electric field gradients in the propaga- 
tion of spin-polarized carriers. 

We note that spin accumulation at the interface be- 
tween regions with different doping concentrations, such 
as seen in Fig.[2l^a), has been previously discussed by Per- 
shin and Privman.— It was argued that injecting spin- 
polarized electrons across a low-to-high doping concen- 
tration interface yields a significant enhancement of the 
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FIG. 4: (color online) Spin density imbalance 5-\y (a), 
spin density polarization P„ (b), and electric field distri- 
bution (c) for a spin and charge inhomogeneous structure 
as depicted in Fig. [T] for different channel lengths, L = 
0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.8, 1.0, 1.8, 3.0 /xm. Curves are offset 
such that a;p = in all cases. Parameters are Nl = 10^^ 
m~^ Nr = 10^^ m-^ Ts/ = 1 ns, = 0.1 ps, Vb = -0.3 V. 
Arrows indicate the position at Xd- The thick, solid line in 
the figures corresponds to the charge homogeneous case. 



spin polarization. However, the increase is seen only 
for (5||, i.e., the spin density imbalance^ which the au- 
thors considered in their work. This is, however, only 
due to the overall total charge density increase as dis- 
cussed above. A closer look and comparison to the total 
charge density, as shown in Fig. [2^ a) showed that is 
suppressed at the interface, which becomes even more 
evident when considering the actual spin density polar- 
ization^ Pn- Moreover, we find that P„ can be enhanced 
for spin-polarized electron injection across a high-to-low 
doping concentration interface, due to the positive elec- 
tric field gradients [see dashed-dotted lines in Fig. [21[c)]. 
We would also like to point out the large difference be- 
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tween the spin density imbalance, and polarization, 
P„, in inliomogeneous systems. Previously in the liter- 
ature, these two quantities have been used interchange- 
ably, which certainly is true for charge homogeneous sys- 
tems. For charge inhomogeneous systems on the other 
hand, depending on the measuring setup at hand, the 
correct quantity has to be studied and discussed for the 
optimization and assesment of spin injection, transport 
and detection schemes. 

We again emphasize that the observed characteristics 
in the polarization, P„, are quite unexpected. Naively, 
one would expect that, e.g., an increase in doping con- 
centration would result in a simple decrease of P„ (which 
is indeed observed), due to the rapid increase in the to- 
tal charge density, n. However, it is not obvious how P„ 
should vary for different doping concentrations, as in the 
setup discussed in our paper, and what the magnitude 
and spatial dependence of P„ should be. Our numerical 
calculations, discussion and analysis above, clearly show 
that the spatial dependence of the spin density polariza- 
tion is quite complicated in the presence of inhomoge- 
neous electric fields. 



significantly change the spin diffusion length. Hence, 
purely electrical effects may have a profound impact on 
spin injection, spin transport, and spin detection, three 
key ingredients for the realization of semiconductor spin- 
tronics. Our findings should be relevant to current exper- 
imental efforts, as well as theoretical studies in the field, 
where we have shown that it is crucial to take into ac- 
count inhomogeneous electric fields self-consistently, in 
particular around doping interfaces as discussed in our 
paper. 
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